A time-dependent variational principle for dissipative dynamics 
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We extend the time-dependent variational principle to the setting of dissipative dynamics. This 
provides a locally optimal (in time) approximation to the dynamics of any Lindblad equation within 
a given variational manifold of mixed states. In contrast to the pure-state setting there is no 
canonical information geometry for mixed states and this leads to a family of possible trajectories 
— one for each information metric. We focus on the case of the operationally motivated family of 
monotone riemannian metrics and show further, that in the particular case where the variational 
manifold is given by the set of fermionic gaussian states all of these possible trajectories coincide. 
We illustrate our results in the case of the Hubbard model subject to spin decoherence. 



I. INTRODUCTION 

One of the main challenges in a quantum mechanical 
experiment is to overcome the interaction of a system 
with its environment. Such interactions lead to decoher- 
ence and often obscure coherent quantum phenomena. 
Recently it has been shown that this vice can be turned 
into a virtue: dissipative processes can be exploited as 
a possible resource for quantum state engineering [T]- 
[3]- Several evolutions leading to non-trivial fixed points 
have now been proposed, including states with non-trivial 
topological properties |4j. Such dissipative engineering 
has opened up a completely new world for us. 

Motivated by the new possibilities offered by dissipa- 
tive engineering there has been renewed interest in under- 
standing dissipative processes in more detail. However, 
this task is complicated by the fact that, as for ground 
states, we can only hope for analytic solutions in very spe- 
cial cases. Therefore, in general, we must take recourse 
to numerical approximation techniques in order to gain 
insight into the physics of a dissipative system. Typically 
the method of choice here is a Monte Carlo sampling al- 
gorithm. Such methods have led to many insights into 
the dissipative systems occurring in quantum optics, but 
have faced limitations when applied to strongly interact- 
ing many particle systems, particularly fermions, due to 
the inevitable sign problem. 

There is, however, another general approach available 
to us, namely the variational method. This method 
has been very successfully applied in the pure-state 
case leading to unparalleled insights into the equilib- 
rium physics of strongly interacting many body systems. 
Further, the elegant time- dependent variational principle 
(TDVP) allows the locally optimal study of nonequi- 
librium unitary dynamics. The power of this method is 
well known in the held of quantum chemistry, where its 
application to the class of Hartree-Fock states is known 
as time- dependent Hartree Fock theory [7J. This tech- 
nique has also been exploited to great effect in the context 
of one-dimensional quantum spin systems in conjunction 
with powerful expressive variational classes such as ma- 
trix product states, a method synonymous with the den- 



sity matrix renormalisation group (DMRG) |S]. 

In contrast to the pure-state case, there is no opera- 
tionally unique way to formulate the variational method 
for mixed states because there is no distinguished mea- 
sure of information distance: there arc infinite families of 
inequivalent distance measures, including examples such 
as the fidelity and the trace distance. This has com- 
plicated the formulation of a mixed-state TDVP, which 
requires knowledge of the geometry of state space. How- 
ever, recent results in our understanding of the infor- 
mation geometry of mixed states allow us to revisit this 
problem. (See, however, [51 [TO] for related variational 
approaches to the von Neumann equation.) 

Thus, in this Article, we formulate the TDVP for 
mixed states in the general case of distance measures 
arising from monotone riemannian metrics. We then 
show that this method, when applied to the variational 
class of fermionic Gaussian states evolving according to 
an arbitrary Markovian CPT (completely positive trace- 
preserving) map p t — £t(Po) are all equivalent to the 
application of Wick's theorem (via Gaussification). Fi- 
nally, we apply this method to the one-dimensional (Id) 
spinful Hubbard model subject to a decoherence process. 

II. A REVIEW OF THE TDVP FOR PURE 
STATES 

In this section we review of the TDVP for pure quan- 
tum states, and explain why this approach cannot be 
immediately applied to the mixed case. First, we present 
the necessary notation. We denote by M. n {G) the set of 
all complex n x n matrices with entries in C. The state 
space of an n-dimensional quantum system is given by 
the set T> n (here viewed as a differentiable manifold) of 
all density operators defined by V n = {p e M. n {^)\p^ — 
Pi P > 0,tr(p) = 1}. In order to formulate a time- 
evolution within this manifold, we have to introduce the 
notion of tangent space. The tangent space T p T> n to T> n 
at any interior point p € T> n can be identified with the 
set {A e M n (G) | At = A,tr(A) = 0} of traceless her- 
mitian matrices (we assume that our processes never in- 
clude any part of the boundary of the manifold). The 



2 




FIG. 1: Time-dependent variational principle with respect to 
a variational manifold of mixed states V = {p(x)|x G R}. A 
mixed quantum state p(x) within a variational manifold V 
evolves in time according to an arbitrary physical process de- 
scribed by a CPT map p t = £t{po)- In general, such a process 
leads out of the tangent space of the variational manifold at 
the point p(x). Hence, we have to "project back" into V using 
an appropriate measure of distance (color online). 



set T> n can be given the structure of a Riemannian man- 
ifold by choosing a positive bilinear form M p (A,B) on 
TpV n for all p G T> n . This supplements us with the no- 
tion of a distance. Throughout wc define a variational 
class to be a submanifold V of T> n parametrized accord- 
ing to V = {p(x) |x G MP}, where the dependence on 
the parameters x J is assumed to be analytic. 

We consider the time evolution of a quantum state p t — 
p(x(t)) G T> n in its most general form, i.e. p t = £t{po), 
where £t is a completely positive trace-preserving (CPT) 
map. Assuming that £ t is differentiable with respect to t 
allows us to write the equation of motion as d t pt — L{p t ), 
where C is the infinitesimal generator of the dynamics. 
For example, £ could describe a Hamiltonian or dissipa- 
tive evolution of a quantum system. An exact solution 
to this evolution is in general hard to find, and so we 
aim at finding an optimal approximation to this evolu- 
tion within the variational class V. 

To this end, we first review the case of Hamil- 
tonian time evolution within the set of pure quan- 
tum states. Here, the variational class of state vec- 
tors is represented as the set {|i/>(x))|x G H D }. 
The time-dependent Schroedinger equation then reads 
& dj\ip(x)) = —iH\i/j(x)), where dj — djdxK Note that, 
in general, the vector H\ip(x)} is not an element of the 
tangent space to V at \ip(x)}, whereas the left side is 
a linear combination of vectors that span the tangent 
space IW X )\V. Thus, in general, there is no exact so- 
lution for x J . The best approximation is given by the 
solution to the minimisation of the information distance 
(here measured using the fidelity) between the left- and 
right-hand sides: min^j W&dj + iH |^(x)) || . The 
minimum can be found by applying an orthogonal pro- 
jection of iH\ip(x)) onto the tangent space, as depicted 
in FigJU 

This discussion immediately reveals why the TDVP 



cannot be directly applied to the mixed-state setting: 
the approximation of the RHS of x^djp = Ct(p) by a 
vector in the tangent space requires a unique notion of 
information distance. In the mixed-state case there is 
no operationally unique answer, since there exist infinite 
families of inequivalent measures, and hence, there is no 
canonical choice of Riemannian metric on T> n . However, 
we explain in the next section a possible solution to this 
problem. 

III. FORMULATION OF THE TDVP FOR 
MIXED STATES 

As we have explained above, there exists no canon- 
ical choice of Riemannian metrics in on the set of 
mixed quantum states. However, it turns out that there 
are several families of Riemannian metrics which nat- 
urally arise from information-theoretic considerations. 
Here the natural condition is that the metric is mono- 
tone, meaning that the norm induced by the bilinear 
form Mi p \ cannot increase under any CPT-map £, i.e. 
M £(p) {£{A),£{A)) < M p (A,A). The reasoning here is 
that the distinguishability of two states infinitesimally 
close to p can never be increased under the action of a 
channel. Remarkably, Petz showed there is a one-to-one 
correspondence between the set of monotone metrics and 
a special class of superoperators Q, p (built in terms of con- 
vex operator functions) : these lead to monotone met- 
rics according to M p (A, B) = (A, fl p (B)) = tr[A^Q p (B)}. 

These monotone metrics now allow an operationally 
motivated formulation of the TDVP for the dissipative 
dynamics generated by dtpt — £{pt) within a given vari- 
ational class V. The setup is identical to the pure-state 
case: we aim to find the optimal trajectory p t G V gener- 
ated by the vector field coming from the optimal element 
A G T Pt V which is closest to £(pt), where we use the 
quadratic form M Pt (A, B) to measure the distance. That 
is, we solve infot eTpt y M Pt (A - £(p t ),A - C(p t )). An in- 
tuitive picture of this last equation is given in Fig. [I] The 
evolution under £ takes us out of the variational mani- 
fold V, and we want to "project back" into V to find the 
state in the variational manifold that is the best approx- 
imation to this evolution. Using the definition of M Pt we 
find that this is equivalent to solving 

inf (A-C{p t ),n pt {A-C{ Pt ))). (1) 

Parametrizing A = v j djpt we can rewrite this infimum 
as in{ veR D v T G p v — v T l p — ljv + c$. The solution is 
given by 

v = g;%, (2) 

(G p ) jk = (d jP (x(t)),n p (d kP (x(t)))}, (3) 

(U- = (a, P (x(t)),n p (/:(p(x(t))))), (4) 

where G p is the pullback metric or Gram matrix. This 
solution gives us the optimal trajectory (locally in time) 
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within V via integration of the equation of motion 

d t pt=v : >{t,pt)d j pt. (5) 

This equation of motion is the first main contribution of 
our Article. Eqs. (2)-(5) can be applied to any varia- 
tional manifold subject to any physical process. In the 
next Section we apply this framework to a concrete ex- 
ample that is relevant in many-body physics. 



IV. THE TDVP FOR FERMIONIC GAUSSIAN 
STATES 

The understanding of fermionic quantum systems is 
of central interest many fields of physics. Fermions are 
building blocks of matter and thus central to some of the 
most fascinating effects known in the theory of many- 
body physics, like superconductivity or the quantum Hall 
effect. However, most problems of interest do not have 
a closed analytic solution, and we have to use numeri- 
cal approximation techniques and we have to use appro- 
priate variational wave functions to obtain insight into 
these systems. In the following we introduce the class 
of fermionic Gaussian states (fGS) that has been suc- 
cessfully applied to solve fermionic many-body problems 
in the pure state setting. We show then how the TDVP 
can be applied using this class of states as our variational 
class Vq and give a numerical example in the last Section. 

In the following we describe fermionic systems in terms 
of N fermionic mode operators aj obeying the canonical 
anti-commutation relations {a£.,a;} = Sf-i- We use the 
equivalent representation in terms of 2N hermitian Ma- 
jorana operators C2j— 1 = a] + <Xj & n d C2j — (— *)( a ] — a j) 
which obey {c kl c{\ = 2S k i. We take as our variational 
manifold the set of fermionic Gaussian states (fGS). 
fGS are those states whose density operator can be ex- 
pressed as an exponential of a quadratic function of 
the Majorana operators, p — Kexp[—^c T Kc], where 
K = —K T G R 2Nx2N . All information about the state 
is encoded in the real and anti-symmetric covariance ma- 
trix (CM), Tki — |tr([cfc, ci]p), due to Wick's theorem: 

i?tr\pc h ...c j2p ] = Pf(T h where 1 < ji < ... < 

dip < 2M and ^j 1 ,...j 2p is the corresponding 2p x 2p sub- 
matrix of r. Pf(rji,...j 2 p) 2 = det(rj li ... jj2p ) is called the 
Pfaffian (see, e.g., [T2] for further details). 

The class of fGS is a natural generalization of the 
variational classes used in Hartree-Fock and BCS-theory. 
Thus, it is combining and extending the most success- 
ful tools in the description of fermionic many-body sys- 
tems and hence allows for a description of a wide class 
of fermionic phases of matter, like superfluids, Mott and 
spin ordered phases. Recently it has also been shown that 
fGS with topological order can be engineered in a cold- 
atom implementation via a local dissipative process [4]. 
Thus, fGS have proven to be a powerful class capable of 
capturing fermionic phases with highly non-trivial prop- 
erties. Every pure fGS is the ground state of a quadratic 



Hamiltonian. Further, fGS remain Gaussian under the 
evolution according to a quadratic Hamiltonian or a dis- 
sipative process with linear Lindblad operators. Using 
this, fGS have allowed the approximation of the ground 
and thermal states of, the time-evolution 13J of, as well 
as the excitation spectra [14] of interacting fermionic sys- 
tems. 

The main ingredient used in all studies exploiting fGS 
is a process known as Gaussification, i.e., the approxi- 
mation of any iV-body correlation function in terms of a 
product of single-particle correlation functions via Wick's 
theorem (see above): 

Definition 1 Let a be a fermionic quantum state. Then 
its Gaussification, pq = Q{u) G Vg> is defined via the 
relation T(pg) — r(er), which is equivalent to the appli- 
cation of Wick's theorem to a. 

In the following, we show that the process of Gaussifica- 
tion is locally optimal in time within the variational class 
of fGS for all a metrics, since every monotone metric can 
be written as a convex combination of them |llj . 

Theorem 1 Let £ t be an arbitrary (differentiate) 
Markovian CPT map defining a time evolution on the 
space of density matrices via p t — £t(po)- Then, the 
optimal approximation of this time evolution within the 
variational manifold of Gaussian states Vg with respect 
to any (convex combination of) monotone metrics of a- 
norm type: 

Q a p (<T) = l(p- a ap a - 1 +p a - 1 ap- a ), (6) 

is obtained via Gaussification. 

Theorem 1 is the second main result of this Article. To 
prove it we obtain a Gaussification of the time evolution 
of p according to the generator C of the CPT map St as 
follows. Let pit + St) = p(t) + St£(p), where St is an in- 
finitesimal time step. The operator p(t + St) — p(t) is not 
necessarily a member of the tangent space T p Vg- How- 
ever, the Gaussified operator Q (p(t+St))— p{t) is, and it is 
therefore a linear combination of the tangent vectors djp, 
i.e. G{p(t+St)) = p(t)+5tJ2 3 vidjp, where G R. Since, 
by definition, Gaussification implies T(Q(p)) — T(p), we 
obtain from the linearity of the map T, the following 
defining condition for the expansion parameters v J : 

tr[c kl c k2 C{p)} =^2v J ti[c kl c k2 djp\. (7) 

3 

Now we show that an application of the TDVP projection 
with respect to any a- norm also leads to Eq. ([7| . For the 
proof of this statement we need the following 

Corollary 1 Let p G Vg be a fGS of N modes, and 
let T p Vq denote the tangent space of p. Then, B = 
{d(j lt j 2 )p = (^p) _1 («c :(1 c : , 2 )}i< : , 1<: , 2 <27v is a hermitian 
basis for the tangent space T p Vg- (Here (ji, j'2) denotes 
a multi index. ) 
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Proof We first claim that if {Q a }a=i is a set of linearly 
independent operators then so is {0" (Qa)}a=i- This fol- 
lows immediately from the fact that Q p is a positive su- 
peroperator so that B is a set of linearly independent 
vectors. 

Next, we show that dim(S) = dim(X , p VG). To this 
end, we determine a basis of T p Vg by applying the most 
general infinitesimal Gaussian transformation on p |12j . 
These are of the from p W pW^ /ti[W pW'], where 
W = e ieZ ^ c ' with Z kl = X kl +iY kl , X kl ,Y kl e R, 
and e <IC 1. This leads directly to the tangent vectors 
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[R) 
kl 



= [p,c k ci] and A k y = i{p,c k ci} 



2T kl p. In 



order to determine the number of linearly independent 
tangent vectors, we work in the basis c k = J2i^kici, 
where OO t = 1, so that p is in its standard form, 



P = rij=i |C- + i^jC2j-iC2j), where Xj € (-1,1) 
(i.e. p has no pure subspace). Then, the tangent 
vectors are readily calculated. For all 1 < k < I < N 



we obtain A. 



(R) 

2fc.22 



-2i(Xic 2k c 2 i- 



X k C 2k -lC 2 l)p k l, 



A 2k,2i-i = 2i(AiC 2 fcC 2 ; 
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X k C2k-lC2l-l)Pkh ^2fc-l,2/-l 

2«(A ; c 2fe _ 1 c 2 i + X k c 2k c 2 i-i)p k i and A^f > _ 12l 
2i(-A/c 2fc _ic 2; _i + \ k C2kC2i)pkh where p ki 
n,^fcj |(1 + *^j%'-i9y)- Further, we find 



^2k,2! = 2i(c 2fc C 2 ; 



that 



XkXlC2k-lC2l-l)pkl, A 



A 



CO 

2k,2l-l 



2«(c 2 feC 2 /_i + AfcA/c 2 fc_ic 2 i)/3fci, ^ 2 fe-i, 2 ;-i 

2i(c 2fe _ic 2; _ 1 - X k Xic 2 kC2i)pki and ^4fc-i,2i = 
2«(c 2fe _ic 2; + X k Xic 2k c 2 i~i)pki- Thus we obtain, for all 
A/c,i G (—1, 1), the four linearly independent basis vectors 
ic 2 k-\c 2 i-iPkh ic 2 k-\C2iPki, ic2kC2i-iPki and ic 2k -C 2l p M . 



If k = I. we obtain A 



(i) 

2k-l,2k 



A. 



(R) 

2fc-1.2fc 



2i(l-Xl)c 2k -ic 2k p k , while 

= 0. Thus the dimension of the tangent space is 
dim(T p V G ) = 4N(N-l)/2 + N = 2N(2N-1) = dim(K). 

Finally, we show that the operators d(j 1 ,j 2 )P are her- 
mitian. This follows immediately from the fact that 
^p((d(n^)PV) = ^p(9(n,h)p), and that 0« is invert- 
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With Corollary 1 in hand we now directly apply the 
TDVP projection in our special basis, following Eq. 
We see that 

( G p)0"U2),(fei,fc 2 ) = ( d {n,n)P^ a p( d (kiM)P)) 
= tr [%ij2)P c fci c fc 2 ], 
( l p)(n,32) = ( d Uu32)P> n p( C (p))) = tr[£(p)c fel c fc2 ], 

and we arrive at Eq. Q. This proves the equivalence of 
Gaussification and the application of the TDVP to the 
variational class of Gaussian states. 



V. EXAMPLE 

In the following we apply our approach to the Id spin- 
ful Hubbard model with repulsive interactions subject to 
a magnetic field. The evolution of the system is described 
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FIG. 2: a) Anti-ferromagnetic order (AF) in the ground state 
po of the Hubbard model for u = 4, p, = —2. The steady state 
p s shows a ferromagnetic order (FM). b) Comparison of the 
real-time evolution (time t in units of 1/k) of the dissipative 
process with the Gaussified version. We present the differ- 
ence of the CM and the two states dF(t) — ||r pG ( t ) — r p ( t )||2, 
dp(t) = ||pg(£) — ) 1 j 2 - c) Purity, p p = tr[p 2 ], for the exact 
(solid) and the Gaussified (dashed) evolution, pa- d) Evolu- 
tion of the particle number n for the two spin states a s =t,i 
for the real (p) and Gaussified (pa) process (color online). 



by a Lindblad equation dtp = — i[H, p] + kJ2 x JxPjl 
\{iiix,p\, where 

H = J ^ a l,o- 3 a x+i,cr 3 + u } j n x ,fn x>i , + 



Jx ^X^^ X :^' 

Starting from the ground state of the Hubbard Hamilto- 
nian we expect that the external noisy magnetic field 
modeled by the operators j x drives the system to a 
completely spin-polarized state. The time scale of this 
process depends on the ratio between the decoherence 
strength k and the parameters of the Hubbard Hamilto- 
nian. 

We consider a system of L = 4 sites with periodic 
boundary conditions at half filling and consider an inter- 
action u — 4, a chemical potential p = —2 and k = 1, 
where we take the hopping J to be the energy scale. The 
unique ground state po has anti-ferromagnetic (AF) order 
(see Fig{2^, (squares)). We implement a real-time evolu- 
tion p(t) of the dissipative process, arriving at a unique 
steady state p s which is completely spin polarized (fer- 
romagnetic order (FM)) (Fig. [2^, (circles)) and given by 
rL=i a LtlO)- I n or der to measure how well the Gaussified 
evolution given by pq (t) approximates the exact dynam- 



ics, we present the deviation dT(t) = \\T pG / t -) — r p( - t )|| 2 of 
the CM of the real and Gaussified process and the dis- 
tance between the two states dp{t) = Wpoit) — p(t)\\2 in 
Fig.[2]D. We find that, as one might expect, on short time 
scales the Gaussified dynamics takes a different path from 
the exact evolution, but quickly coincides with the exact 
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evolution for intermediate and long time scales. This 
can be explained by the fact that the ground state of the 
Hubbard model is not a Gaussian state, so that the ini- 
tial state for the Gaussified evolution is a large distance 
from the ground state. 

In order to obtain more insight into the Gaussified evo- 
lution we compare the time dependence of some physical 
quantities with the exact evolution. In Fig.[2]c we present 
the purity p p = ti[p 2 ] for the exact evolution (solid) and 
the Gaussified process (dashed). We see that we end up 
in a pure state with only spin up particles (Fig. [2ji) in 
the limit t — > oo in both cases. 



VI. SUMMARY 

In summary, we have extended the time-dependent 
variational principle to the mixed-state setting, provid- 
ing locally optimal (in time) equations that allow for an 



approximation of any Lindblad dynamics, given a vari- 
ational manifold of mixed states and some information 
metric. In the case of fermionic Gaussian states we have 
proven that all a-metrics lead to the same dynamics in 
the space of density matrices, which can equivalently be 
obtained via an application of Wick's theorem at each 
time step (Gaussification). Thus, this method can easily 
be applied to large systems in any dimension and geome- 
try, providing a powerful numerical tool for a variational 
study of dissipative dynamics for mixed quantum states. 
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